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METHODS FOR PROCESSING 
ELECTROCARDIAC SIGNALS HAVING SUPERIMPOSED COMPLEXES 



This patent application claims the benefit of priority under 35 U.S.C. § 1 19 of 
1 5© U.S. Provisional Application Serial No. 60/247,269, filed November 10, 2000, entitled "Method 

O for Viewing and Comparing ECG Signals Having Superimposed Complexes," and U.S. 

in 

T. Provisional Application Serial No. 60/295,217, filed June 1, 2001, entitled "An Algorithm to 

Q 

Measure T-wave Subtraction Quality/' the entirety of both of these applications being hereby 
U incorporated by reference. 
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Field of the Invention 

This invention relates to a method for processing electrical signals obtained from 
the heart and, more particularly, to a method for processing electrocardiac signals having 
superimposed sub-component complexes to enable tracking of native, paced, and derived beat 
25 signals. 
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Background of the Invention 

Certain cardiac arrhythmias are triggered or initiated from a site in the heart tissue 
other than the sinus node. These arrhythmias are generally classified as being "focal" in nature. 
Treatment of focal arrhythmias generally involves locating the arrhythmogenic site and ablating 
it. One method for regionally locating the focal site is the use of a diagnostic 12 Lead ECG. The 
12 Lead can be used in conjunction with pacing via a roving intracardiac catheter to pace map 
the heart. The theoretical basis of this method assumes that the paced 12 lead ECG will appear 
identical to the non-paced ECG if the cycle length (i.e., paced heart rate) and pacing site matches 
the non-paced heart rate and focal site of origin. 

One problem with this method (in current practice) is the subjectivity involved in 
visually comparing a non-paced 12 Lead ECG to a paced 12 Lead ECG. 

A second problem is the time consuming nature of the procedure in which, 
typically, a spontaneous ectopic beat is recorded and printed on paper. A roving mapping 
catheter is positioned at a likely site of ectopy, pacing is initiated, a recording is made, a printout 
is generated and a visual comparison is made by aligning the printouts from the spontaneous and 
paced beats over one another. This process is repeated in an iterative manner until the physician 
determines that a good match between the spontaneous ectopic beat and the paced beat is found. 

A third problem arises when multiple arrhythmogenic foci are present and each 
focus produces a variant on the 12 Lead ECG. Better discrimination between these foci would 
be advantageous during pace mapping as well as during other EP procedures. (Ref - Throne RD, 
Jenkins JM, Winston SA, et al. "Use of tachycardia templates for recognition of recurrent 
monomorphic VT " Comp. Cardiology 1989:171-174.) 



A fourth problem involves the superimposition of the P-wave and T-wave 
components of the ECG. The electrocardiogram typically includes an initial impulse, termed the 
P-wave, emanating from the atria, followed by what is termed the QRS complex, emanating 
from the ventricles, which is followed by a T-wave resulting from repolarization of the ventricles 
(Fig. 1). Thus, a heart beat begins with the P-wave and ends with the T-wave, and the next heart 
beat begins with another P-wave. 

The P-wave can be a valuable tool used by clinicians to diagnose the condition of 
the heart. Thus, clinicians will often monitor an electrocardiogram (ECG) of the heart to aid in 
the diagnosis of atrial and ventricular arrhythmias. This can be done in various ways, such as by 
monitoring the 12 Lead (surface) ECG in conjunction with observing the bioelectric activity 
recorded on intracardiac electrodes carried by a transthoracic catheter. 

In some focal arrhythmias the atrial heart tissue begins to beat very rapidly as the 
focal origin moves from the sinus node to an ectopic site. Sometimes this higher heart rate is 
sustained over three or more beats and is termed a tachycardia. Other times the higher rate is 
intermittent and may be as short as one heart beat. In either case, the first beat of the atrial 
arrhythmia is usually initiated by what is termed a Premature Atrial Contraction ("PAC") which 
can result in the P-wave of a successive heart beat overlapping with the T-wave of the preceding 
beat (Fig, 2). Not only is this a physiologically compromised state for the heart to be in, but the 
clinician can no longer use the P-wave to diagnose the heart because it is obscured by the T- 
wave. 

Accordingly, it will be apparent that there continues to be a need for a method that 
allows a clinician to Pace Map more effectively and in addition monitor the P-wave of a patient's 



heart beat, even when the P-wave is overlapping with a preceding T-wave. The instant invention 
addresses these needs. 

And while T-wave subtraction is a useful method in electrophysiology procedures 
to unmask the ECG P-wave morphology of a PAC by subtracting a QRS JT template from a 
5 PAC, ECG baseline drift caused by respiration or body movement may cause certain variations 
on the results of T-wave subtraction. Thus, a further need remains in the art to quantitatively 
measure the quality of T-wave subtraction results, among other reasons to monitor the respiration 
variations on T-wave subtraction. The instant invention addresses this need as well. 

i ■, 

log Summary and Objects of the Invention 

.0 

tfl The present invention, in certain aspects, provides a medical practitioner with a 

N| computerized method for objectively and efficiently performing real time pace mapping and 

Q 

: x other cardiac analyses, through the processing of incoming electrical signals which represent 

fl heart activity to display a derived P-wave without any overlap with a preceding T-wave during a 

3 

15 q PAC, and to allow the practitioner to objectively compare derived P-waves to determine if they 
are emanating from the same focus. As a direct consequence of the cardiac signal processing of 
the present invention, otherwise masked signals and correlations are identified among heart beats 
and segments of heart beats through calculations on acquired signals and/or derivations of new 
signals. The practitioner can be guided through visual aids such as bar graphs and overlaid 

20 cardiac signals of the quality of signal matches. These signal matches can assist in diagnosing a 
patient and in the effectiveness of an ongoing treatment, for example, an ablation procedure. 
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f Due to timing and amplitude relationships among beats of a heart, there is the 

possibility that individual waveforms can be obscured or hidden. If a singular, unadulterated sub- 
component waveform is identified, and if this sub-component has similar timing characteristics 
that allow it to be synchronized with the composite waveform, then a subtraction process can be 
5 performed in accordance with an aspect of the invention to thereby derive the other sub- 
component waveform(s). Sub-component waveforms, either derived, native state, or pace 
induced, can be quantitatively compared to one another using correlation analysis. This analysis 
may be done retrospectively or in real time, 
y. More specifically, the present invention provides systems, programmed machines, 

1 0 q and methods that permit superior signal processing over prior art electrophysiology signal 
IH processors and can achieve this using a standard 12 lead ECG.. 

2 In accordance with one aspect of the invention, a system for tracking ectopic beats 

J s comprises a signal sensing unit, a signal processor, and an output device. The signal sensing unit 
p*! is configured to capture a first ECG signal. The signal processor is connected to receive the first 
15 q ECG signal from the signal sensing unit and is configured to permit a user to mark a begin point 

Li 
teas*. 

and an end point of the first ECG signal for use in defining a waveform segment as a reference 
template, to acquire data from multiple leads, and to identify a best fit between the reference 
template and the acquired data using a correlation coefficient calculation. The output device 
presents the identified best fit. 
20 In accordance with another aspect of the invention, a system for deriving a p- 

wave signal from a premature atrial contraction ("PAC") beat comprises a signal sensing unit, a 
signal processor, and an output device. The signal processor is connected so as to receive 
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„ electrocardiac signals from the signal sensing unit and is configured to process the electrocardiac 

signals so as to derive the p-wave signal from the PAC beat. The output device presents the 
derived p-wave signal 

In a particular embodiment of the foregoing system, the processor is configured to 
5 execute the steps of: (a) selecting a QRS-T segment of a reference ECG signal; (b) permitting a 
user to mark a begin point and an end point of the selected ECG signal; (c) defining a reference 
template as being a waveform segment between the marked begin and end points of the selected 
ECG signal; (d) acquiring the PAC beat at the signal processing unit from multiple leads 
i (preferably with no more than 12 leads); and (e) processing the PAC beat so as to derive the p- 

10pi wave signal 

Q 

[| In accordance with still another aspect of the invention an electrophysiology 

SI computer system includes a processor that is configured to derive a p-wave signal hidden within 
a a premature atrial contraction ("PAC") beat. The processor executes the steps of: (a) selecting a 
QRS-T segment of a reference ECG signal; (b) permitting a user to mark a begin point and an 
15 f 2 end point of the selected segment of the reference ECG signal; (c) defining a reference template 
^ as being a waveform segment between the marked begin and end points of the selected segment 
of the reference ECG signal; (d) acquiring the PAC beat at the signal processing unit from 
multiple ECG leads; and (e) processing the PAC beat so as to derive the p-wave signal 

In a particular embodiment of the foregoing system, the processor utilizes a 
20 correlation coefficient calculation to effect a subtraction of the reference template from a 

predetermined segment of the PAC beat. In more particular embodiments, the processor is 
configured to compare derived p-waves from multiple beats to one another, to indicate or infer a 



common focal origin among several derived p-waves, to predict the most likely site of the origin 
of a focus using a (preferably 12 lead) library of p~waves of known focal origin, to derive paced 
p-waves for comparison to spontaneous p-waves, to determine an integral value of the QRS area 
of a derived p-wave signal, to normalize any integral values over a length of the derived p-wave 
signal, to process the QRS segment of a beat separately to arrive at further determinations 
concerning the heart beat data, and to perform combinations of the foregoing. 

In accordance with still another aspect of the invention an electrophysiology computer 
system includes a processor that is configured to execute steps substantially in the same manner 
as the processor that derives a p-wave from a PAC beat, but more generally is configured to 
derive a non-synchronous subcomponent from a first heartbeat signal having a composite 
waveform which includes a synchronous subcomponent overlapping the non-synchronous 
subcomponent. The processor executes the steps of selecting a synchronous subcomponent of a 
second heartbeat signal which corresponds to the synchronous subcomponent of the first 
heartbeat signal; permitting a user to mark a begin point and an end point of the selected 
synchronous subcomponent; defining a reference template as being a waveform segment 
between the marked begin and end points of the selected synchronous subcomponent; acquiring 
the composite waveform of the first heartbeat signal at the signal processing unit from multiple 
leads; and processing the composite waveform beat so as to derive the non-synchronous 
subcomponent. 

In accordance with yet further aspects of the invention, a method for tracking 
ectopic beats through template matching is described which includes the steps of: (a) capturing a 
first ECG signal in a signal processing unit; (b) permitting a user to mark a begin point and an 



end point of the captured first ECG signal; (c) defining a reference template as being a waveform 
segment between the marked begin and end points of the first ECG signal; (d) acquiring data at 
the signal processing unit; and (e) using a correlation coefficient calculation on the acquired data 
to identify a best fit between the reference template and the acquired data. 

In accordance with further aspects of the invention, a method for deriving a p- 
wave signal from a premature atrial contraction ("PAC") beat is described which can assist a 
person in diagnosing a heart. This method includes the steps of: (a)selecting a QRS-T segment 
of a reference ECG signal; (b) permitting a user to mark a begin point and an end point of the 
selected segment of the reference ECG signal; (c) defining a reference template as being a 
waveform segment between the marked begin and end points of the selected segment of the 
reference ECG signal; (d) acquiring the PAC beat at the signal processing unit from multiple 
leads; and (e) processing the PAC beat so as to derive the p-wave signal. 

In a particular embodiment of the foregoing methods, the PAC beat is processed 
using a correlation coefficient calculation to effect a subtraction of the reference template from a 
predetermined segment of the PAC beat. Also, the foregoing methods can include the additional 
steps of: comparing derived p-waves from multiple beats to one another; indicating or inferring a 
common focal origin among several derived p-waves; predicting the most likely site of the origin 
of a focus using a (preferably 12 lead) library of p-waves of known focal origin; deriving paced 
p-waves for comparison to spontaneous p-waves; determining an integral value of the QRS area 
of a derived p-wave signal; normalizing any integral values over a length of the derived p-wave 
signal; processing QRS segment of a beat separately to arrive at further determinations 
concerning the heart beat data, and performing combinations of the foregoing steps. 
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Further methods according to still further aspects of the invention include the 
determination of integrals concerning a section of the QRS_T segment and the processing of 
those integrals. A QRS segment integral can be used as a measure of the QRS residue, which is 
an indicator of the alignment or synchronization quality between the template QRS and the PAC 
QRS Furthermore, baseline drift can be monitored as a change of the QRS absolute peak 
(integral) value percentage between the template and the PAC. These methods are implemented 
by suitably configured computer processors. 

Yet a further method in accordance with another aspect of the invention proceeds 
in substantially the same manner as when deriving a p-wave from a PAC beat, but more 



10X| generally includes the selecting the synchronous subcomponent of the heartbeat signal, 

Q 

jjj permitting a user to mark a begin point and an end point of the selected synchronous 

subcomponent, defining a reference template as being a waveform segment between the marked 
begin and end points of the selected synchronous subcomponent, acquiring the composite 
p, waveform at the signal processing unit from multiple leads, and processing the composite 



15 m waveform beat so as to derive the non-synchronous subcomponent. 

In accordance with a further aspect of the present invention, a template 
optimization method is disclosed which dynamically employs different templates. QRS beats 
that precede or follow a PAC can be selected manually or by action of a programmed machine in 
selecting and setting a new template for use in subsequent calculations. The method is 

20 implemented by suitably configured computer processors. 
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Other aspects, features and advantages of the invention can be more clearly 
understood from the following detailed description of exemplary embodiments and 
accompanying Drawing Figures. 



Description of the Drawing Figures 

Fig. 1 is a schematic diagram of a normal heart beat; 

Fig. 2 is a schematic diagram of a pre-mature atrial contraction (PAC); 

Fig, 3 is a schematic diagram of the T-Wave subtraction; 

Fig. 4 is a block diagram of a system programmed to practice a method in 
accordance with a preferred embodiment of the present invention; 

Fig 5 is a flow diagram showing the process for template matching in accordance 
with the preferred embodiment. 

Fig. 6 is a flow diagram showing the process for T-wave Subtraction in 
accordance with the preferred embodiment. 

Fig. 7 is a representative computer display interface for T-wave subtraction that 
can be displayed to an operator. 

Fig. 8 is a representative computer display interface for template matching that 
can be displayed to an operator. 

Fig. 9 illustrates a methodology for determining the integrals of a section of the 
QRST segment after the subtraction process. 

Detailed Description of the Preferred Embodiments 

10 



To promote an understanding of the methods that can be practiced in accordance with 
preferred embodiments of the present invention, several pertinent aspects are discussed below 
under respective headings. 

Template Matching/Pace Mapping 

Any recorded ECG waveform can be used as a reference to compare to another 
recorded ECG waveform or to a real time ECG waveform. The comparison is performed in a 
two step process in which first a reference template is selected by the user to describe the 
beginning and end of an ECG waveform segment to be used as a comparison template. Next the 
user selects the region of data to be used for comparison-either from pre-recorded data or from 
the real time data stream. A suitably configured computer processor can find the best match 
against the reference template over the region specified, or in the case of real time analysis, find 
the best match updated over a defined period of time, for example every second. The criteria for 
"best match" utilizes a correlation coefficient calculation across all twelve leads of the ECG and 
finds the best alignment. A visual display showing the aligned reference beat (template) overlaid 
on the beat undergoing analysis give the user feedback as to the closeness of the match. A 
correlation coefficient calculated for each ECG lead gives a quantitative indicator of the match. 
A composite average is also calculated and is displayed in a unique color enhanced bar graph 
indicator which is especially useful when real time template matching is being performed. The 
composite average can be updated as a moving average over a preselected number of beats. 

Template matching may be used to compare two spontaneous beats or it can be 
used to pace map, i.e., to compare a paced beat to a spontaneous beat. A Region of Interest 
(ROI) indictor can be manipulated by the user to exclude certain portions of the waveform from 
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analysis. This is useful during pace mapping where pacing artifacts on the surface leads can be 
excluded from the region of analysis. The ROI indicator can also be used to specify a preference 
for T-wave or P-wave matching as they are oftentimes morphologically very similar. 
T-Wave Subtraction 

5 In one embodiment of the present invention, a method is provided whereby an 

ECG having an overlapping P and T wave is processed to remove the T-wave and thereby 
display the P-wave without any overlap, so that a clinician may observe the P-wave when 
performing a diagnosis of the heart, 
j . See Figure 1 which describes a normal ECG over three beats in which distinctive 

WS52 

O 

1 0 15 P and T-waves can be identified. Figure 2 shows a rhythm in which the P- wave from the third 



beat (P') arrives early and is obscured by the T-wave from the second beat. This results in what 
is termed a P on T complex, and is referred to as a QRS-T-P' in the figure. 

In general, according to the method, the QRS-T segment of a beat that lacks a 
PAC is selected as a template. This template is subtracted from the QRS-T-P' signal in the PAC 



15 p to be studied yielding the P-wave. The QRS-T signal used as the template may be from a single 
beat or it may be derived from an average of multiple beats. The QRS-T signal (or average) used 
as the template is selected so that the preceding QRS-QRS interval is equal (or nearly equal) to 
the QRS-QRS interval immediately preceding the QRS-T-P' signal to be studied. Preferably, the 
beat immediately preceding the PAC can be used for the selected QRS-T template as the cycle 
20 length and hemodynamic conditions of this beat are the closest to those of the succeeding beat 
that contains the PAC and P on T complex. (See Figures 2 and 3.) 



12 



The QRS complex is used as a means to synchronize and align the QRS-T 
template and the PAC beat for subtraction. The alignment is automated by the algorithm for the 
best match based on the composite correlation coefficient across the 12 Lead ECG. The 
practitioner has the option of shifting the template match left or right on a sample by sample 
basis with the resulting composite correlation coefficient updated at each new position. The 
practitioner also has the option of choosing the previous or following QRS-T segment as the 
reference template. The software will automatically locate the previous or following beat based 
on the current reference template and use the corresponding QRS-T segment of that beat as the 
new reference template in the calculation of derived P-waves. 

Different display views showing the derived P-wave, alone, or overlaid with the 
original PAC beat or reference template are available as an aid to the practitioner. 

P-waves that have been derived using the T-wave subtraction method can be 
signal processed further to remove unwanted artifacts caused by respiration or noise. 

3. Template Matching of Derived P- Waves 

Once one has a derived P-wave identified from the tachycardia or premature atrial 
beat (PAC), one can compare this derived P-wave with a previously captured reference template. 

3a. More specifically, one or more spontaneous P-waves may be identified 
using the subtraction method described above and compared with one another using a correlation 
waveform analysis. This can be used to determine if the spontaneous P-waves have the same 
focal origin. This can be done in real time or in review from recorded data. 

3b. In addition, one or more derived spontaneous P-waves may be identified and 
compared to a library of P-waves of known focal origin to predict the most likely site of origin. 
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3c. In addition, once a derived spontaneous P-wave is identified by the T-wave 
Subtraction method as described above then the practitioner can begin atrial pace mapping 
following the Template Matching/Pace Mapping method also described above. The roving pace 
mapping catheter is maneuvered within the atria (or adjacent vessels such as the pulmonary 
veins) until the derived paced P-wave is nearly identical to the derived spontaneous P-wave. 
This comparison of derived P-waves may be done on pre-recorded data or in real time. 

More generally, two or more waveforms X,Y, .., may form a composite waveform 
that due to timing and amplitude relationships causes the individual waveforms to be obscured or 
hidden. The composite waveform includes includes a synchronous subcomponent overlapping a 
non-synchronous subcomponent. If a singular, unadulterated sub-component waveform (e.g. X 
or Y) can be identified, and if it has similar timing characteristics that allow it to be synchronized 
with the composite waveform (i.e., this identified subcomponent is the synchronous 
subcompent), then it can be subtracted from the composite waveform to derive the other sub- 
component waveform(s) (i.e., the non-synchronous subcomponent(s)). Sub-component 
waveforms, either derived, native state, or pace induced, can be quantitatively compared to one 
another using correlation analysis. This analysis may be done retrospectively or in real time. 
One of skill in the art will appreciate that a number of algorithms can be used to compare 
waveform shape, including, but not limited to bin area methods and integrals; any of these 
methods can assist in the goals of aligning synchronous components of composite waveforms 
and/or comparing the derived results. 

A method in accordance with this more general teaching proceeds generally as 
outlined above. Specifically, this method proceeds in substantially the same manner as when 
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deriving a p-wave from a PAC beat, but more generally includes the selecting the synchronous 
subcomponent of the heartbeat signal, permitting a user to mark a begin point and an end point of 
the selected synchronous subcomponent, defining a reference template as being a waveform 
segment between the marked begin and end points of the selected synchronous subcomponent, 
acquiring the composite waveform at the signal processing unit from multiple leads, and 
processing the composite waveform beat so as to derive the non-synchronous subcomponent. 

Referring now to the drawings, and particularly to FIG. 4, there is shown a system 
10 for receiving and processing electrical signals according to one illustrative embodiment of the 
present invention. In one illustrative embodiment, the system 10 includes a signal sensing unit 
12, which may take different forms, such as a standard 12 lead ECG, intracardiac lead, or 
combination thereof. The signal sensing unit is electrically connected to a signal processing 
device 14, which receives the sensed signals from the unit 12 and processes the signals, as is 
described in more detail below. The signal processing device ("signal processor" or "processor") 
14 is preferably connected to a suitable display 16, which will present the processed signals to a 
clinician or other interested person. Information can be stored and recalled from a storage device 
18. Preferably the signal processing device 14 and display 16 comprise the EP LabSystem 
(trademark) of C.R. Bard, Inc., Murray Hill, New Jersey, or the like. The EP LabSystem 
(trademark) supports a variety of data gathering and processing functions that are standard in 
electrophysiology procedures, and can have its hardware (namely, processor 14) configured to 
implement the subtraction and derivation methods set forth above, for example, through software 
(e.g., modules, procedures, functions, or objects) or firmware. The processor 14 communicates 
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with a memory or storage 18 which configures the processor to implement the subtraction and 
derivation methods above (as well as the integral techniques described below). 

In one illustrative embodiment, the special features of the system of the present 
invention are implemented, in part, by a processor using program information stored in a 
memory of the signal processing device 14. The processor 14 can access one or more files, as 
necessary, to implement the required functions, as described in greater detail in connection with 
FIG. 5 and FIG 6. 

Referring now to FIG. 5, the operation of the signal processing device 14 of the 
present invention is described in conjunction with the above structural description of the system 
lOg 10. As illustrated in FIG. 5, the process begins when a clinician desires to create a reference 
P template, and this occurs by capturing a reference ECG signal, as indicated at step 502. 
X, Preferably, the reference ECG signal is captured using a standard 12 lead device and/or one or 

more intracardiac leads. As explained above in connection with Fig. 2, the QRS-T signal 
hk components of a beat which does not exhibit P- on T-wave are selected as a template and it is 
15 S|- this set of electrocardiac signal components that is captured at step 502. Such a beat can be 

captured in sinus rhythm or during a focal arrhythmia such as a tachycardia. Furthermore, it is 
contemplated that the reference template results from signals captured either at the surface, from 
intracardiac leads that can be placed in a variety of locations within the heart, or a combination 
of signals from surface and intracardiac leads. The QRS-T signal that is used as the template can 
20 be captured from a single heartbeat or may be a signal derived from an average of multiple heart 
beats. 



16 



At step 504, beginning and end points of the reference template are marked by the 



clinician using an interface to the signal processing unit 14. The marked points define the 



segment of the ECG waveform to be used as a comparison template. 



At step 506, the clinician selects whether recorded or real-time data is to be used 



5 in the template matching analysis, (This step can be performed at any time prior to the 

waveform matching analysis at step 508, for example, prior to performing steps 502 and 504.) If 
recorded data is to be used in the template matching analysis, then a specified region of pre- 
recorded data is provided to the signal processing unit for comparison to the reference template. 
On the other hand, if real-time data is to be used in the template matching analysis, a stream of 
10pj data from ECG leads is provided to the signal processing unit 14 over a defined period of time 
for comparison to the reference template. 



15 H persons) the result of the template match. The results can be shown qualitatively as 

^* superimposed ECG waveform signals, namely, the reference beat (template) overlaid upon the 
beat under analysis to show the degree of alignment therebetween, or quantitatively as a 
correlation coefficient calculated for each ECG lead. Preferably, a composite average is also 
calculated and displayed. This is illustrated in the computer display shown in FIG 8. 
20 At step 512, a test is made to determine whether the user had selected real-time 

processing at step 506. If so, then the flow loops back to step 508 to again perform the template 
matching analysis and to update the display accordingly. Otherwise, if previously recorded 

17 




alignment between the selected region or time period and the reference template. 



At step 510, the display 16 is updated to indicate to the clinician (or other 



At step 508, the signal processor 14 finds a "best match," in other words, a best 



segments are being analyzed, the user is given the option to save the analysis (as tested at step 
514), and the correlation analysis is saved, as indicated at step 516. Real-time analyses can also 
be saved if desired. 

Referring now to FIG. 6, the operation of the signal processing device 14 of the 
present invention is described in conjunction with the above structural description of the system 
10. As illustrated in FIG. 6, the process begins at step 602 when a clinician captures a PAC and 
desires to subtract a QRS-T reference template from the PAC. The QRS-T reference template is 
marked by the clinician at step 604 (as described above) and a region encompassing the PAC is 
selected by the clinician at step 606 for analysis. The QRS portion of the reference template is 
aligned for best fit with the QRS complex immediately preceding the PAC at step 608. When 
the best fit is found, the processor 14 subtracts the QRS-T reference template from the QRS-T-P' 
segment of the PAC at step 610. 

The difference is the derived P-wave which is output to the display 16 at step 612. 
This is illustrated in the computer display shown in FIG 7, in which the leftmost window 
displays the selected QRS-T reference template between two vertical lines (one dashed line prior 
to the 14 second mark at the top (highlighted by an arrow), and a second solid line just after the 
14 second mark). The rightmost window shows the original PAC waveform with the derived P- 
wave overlaid on top of the portion of the ECG which occurs in the first 15 seconds. The 
overlaid and derived P-wave appears as a second graph superimposed over the ECG signals. 
Visual aids can be provided to automatically align and overlay waveforms for visual comparison 
on a computer display or a printout. 
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Figure 8 illustrates an exemplary display for template matching (without 
subtraction) that can be displayed to an operator. The leftmost window displays markers which 
signify the presence and use of the reference template; the reference template beginning at the 
leftmost vertical line (highlighted by the arrow) and ends at the second vertical line. In this 
example, the reference template marks the start and finish of a P-wave; however, any waveform 
segment can be used if the region of interest has been marked for use as a template. The larger 
display window to the right shows the correlation value for each channel of the 12 Lead ECG as 
compared to the reference template. The bar graph at the far right is inactive in this example 
because the analysis region is taken from recorded data rather than real-time data gathered during 
a medical procedure. 

The data can be saved, printed or both, if desired, in response to a user input to do 
so, as tested at step 614 and implemented at step 616. 

From the foregoing, it will be apparent to those skilled in the art that the present 
invention provides a method for reliably and efficiently recovering a P-wave from a waveform 
that has overlapping P- and T-waves. Furthermore, the template matching capabilities of the 
invention provide the added benefit of quickly and objectively comparing ECG waveform 
components, in their native or derived state. It should also be understood that the correlation, 
subtraction and derivation methods described herein apply to data that can be acquired from 
conventional 12 lead surface ECG signals as well as intracardiac signals or combinations of both 
surface and intracardiac signals. 

Two waveforms can have a high correlation to each other but still be poorly 
matched in absolute terms due to amplitude variation and drift caused by the effects of 
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respiration. This can be a problem when two waveforms are aligned and then subtracted, one 
from the other. It is for this reason that immediately adjacent beats are usually desirable as the 
reference (QRS-T) and PAC (QRS-T-P'). This is not always possible and is not practical when 
performing real time pace mapping. 

A methodology for monitoring the quality of the T-Wave subtraction is now 
described with reference to Fig. 9. At step 902, a subtraction process (as illustrated in Figs. 3 
and 6 and described above) is performed to subtract a QRS-T template from a PAC (QRS-T-P') 
and thereby derive a waveform. The method of Fig. 9 proceeds by then providing integral 
calculations that enable a number of measurements of interest to practitioners, including, but not 
limited to: measures of QRS residue and the quality of the T-wave subtraction process; measures 
of the baseline drift, if any; and optimization of the selection of templates to be used in the 
subtraction process. 

At step 904, the area of a derived waveform is measured. At step 906, the integral 
value is divided by the length of the derived waveform to normalize its value. In addition, at step 
908, the amplitude of the normalized integral value is measured and displayed as a voltage at the 
ECG channel's input. This voltage value is termed the QRS residue. 

As described earlier, correlation analysis is used to align the QRS segment of a reference ECG 
template with the QRS segment of a PAC beat. Thus a further improvement may use the 
correlation coefficient in conjunction with the so-called QRS residue of the derived waveform to 
give an indication of the quality if the match between two beats chosen for subtraction. 
Together, they provide an indicator of the alignment or synchronization quality between the 
template QRS and the PAC QRS. For a perfect alignment and good subtraction results, the 
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derived QRS segment should be flat indicating a high correlation to the template and the QRS 
residue should be very small indicating a small difference in absolute amplitudes (including 
drift). 

Having thus described preferred embodiments of the present invention, it is to be 
understood that the above described arrangement and system is merely illustrative of the 
principles of the present invention, and that other arrangements and systems may be devised by 
those skilled in the art without departing from the spirit and scope of the invention as claimed 
below. 
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